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Strongly interacting systems are characterized by heavily dressed entities with internal degrees 
of freedom, which, on a local level, can be described in terms of coherent quantum states. We 
examine the modification of these local coherent quantum states when such entities condense into 
a macroscopic coherent quantum state, such as superfluidity. As an example, we consider a sys- 
tem of electrons coupled to local lattice deformations. Significant changes in the phonon clouds 
surrounding the charge carriers occur when the system develops into a spatially phase-locked state. 
The question of localized self-trapped charge carriers (bipolarons) in the normal state becoming 
delocalized upon entering the superconducting phase is discussed in terms of squeezing of the local 
coherent phonon states. Suggestions for experimental verifications of these features associated with 
the lattice dynamics are made. 

PACS numbers: 71.35.Lk, 61.10.Ht, 73.20.Jc 



INTRODUCTION 

Particles which locally strongly interact with their 
environment polarize this latter and finish up in self- 
trapped states. On a local level, the environment sur- 
rounding the particles can be described in terms of a 
coherent superposition (in form of Glauber states) of the 
elementary bosonic excitations of the bare environment. 
As a result, new composite entities form, given by the 
charge carriers and their surrounding clouds of bosonic 
excitations of the environment. It is at present a mat- 
ter of dispute whether in the normal state such compos- 
ite entities can exist as itinerant Bloch-like states or are 
purely diffusive. The retarded interaction between the 
charge carriers and the bosonic excitations of the envi- 
ronment favors the picture of diffusive motion caused by 
a dephasing of the dynamics of the two constituents. 

The question we want to address here is whether this 
process of dephasing can be blocked, leading to a delo- 
calization when such self-trapped composite entities con- 
dense into a macroscopic coherent quantum state. We 
investigate such a possibility in terms of a model of very 
general form, given by bosonic charge carriers strongly 
coupled to local bosonic excitations of the environment. 
The bosonic charge carriers can, amongst others, rep- 
resent electron pairs induced by strong electron-phonon 
or electron-spin fluctuation interaction, electron pairs in 
double valence fluctuation systems (stabilized by their 
ligand environment), Frenkcl excitons (stabilized by their 
interaction with the surrounding lattice deformations) 
and, in a more remote sense, He'* atoms in a bath of 
He^ with whose excitations they interact. Such bosonic 
entities can in principle condense into a superfluid state. 
Our interest here is to show how in such a macroscopic 



coherent quantum state the internal structure of the local 
coherent quantum state of the bosonic excitations of the 
environment is modified. 

We shall investigate this problem on the basis of a 
model of itinerant electrons, coupled to some purely lo- 
cal vibrational modes. For moderately strong and strong 
coupling, it results in self-trapped composite entities, 
given by on-site electron pairs, surrounded by local lattice 
deformations which are described by clouds of phonons 
in local coherent quantum states. In the intermediate 
coupling limit these localized entities can be considered 
as two-particle resonant states of the underlying system 
of itinerant electrons. They form out of pairs of such un- 
correlated itinerant electrons and ultimately decay into 
them. 

THE BOSON-FERMION MODEL 

A model Hamiltonian which captures such physics is 
given by the boson-fermion model (BFM) 

^^BFM = {D - n.i^ - t ^ c+ 

+{Ab - 2^i) J2 yt + \\ +vY.[ptc,^c,^^ + pic+c+;\ 

-huJoaY^ (pi + i j (a, + a+) + hwoY^ U+a, + ^ j . 

(1) 

Here denote the creation and annihilation operators 
of hard-core bosons (characterizing the self-trapped local- 
ized electron-pairs) on some effective sites i, which should 
be understood as of being made up of small clusters of 
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atoms rather than single atomic sites. Such hard-core 
bosons have spin-i commutation relations [Pi^,pt]_ = 

2pl and = ^ij- c|j^^ denote the fermionic oper- 

ators referring to annihilation (creation) of itinerant elec- 
trons with spin u, ni„ = c^c^o- being the number opera- 
tor for such fermions. af'^ denote annihilation (creation) 
operators of the excitations of the environment, related to 
the local lattice displacements Xi = {ai + al) / ^j2MujQ/h 
where M is some atomic mass characterizing the effective 
sites, and wo is the frequency of the local deformations of 
the lattice. The bare hopping integral for the itinerant 
electrons is given by t, the bare electron half bandwidth 
hy D — zt (which will serve as energy unit), z denoting 
the coordination number, a denotes the coupling of the 
hard-core bosons to the deformations of the surrounding 
medium and v the exchange coupling between the bosons 
and the pairs of itinerant fermions. The bare boson en- 
ergy level is denoted by A^. The chemical potential ji 
is common to fermions and bosons such as to guarantee 
overall charge conservation. 

The model defined by the Hamiltonian (1) was in- 
troduced shortly after the proposition of the concept of 
bipolaronic superconductivity (appropriate to the ex- 
treme strong coupling anti-adiabatic regime) in an at- 
tempt to extend this concept to the regime of intermedi- 
ate electron-lattice deformation coupling. This model is 
based on a conjecture, rather than on a derivation from 
an underlying basic electron-phonon Hamiltonian, such 
as for instance the Holstein model 

We shall in the following develop the ideas which had 
initially led us to this conjecture. Let us start from a sim- 
ple Holstein model with purely local Einstein modes for 
the lattice vibrations. The corresponding Hamiltonian is 
given by 

H = {D- n)'^n,„~tY^ (c+ c^^ + h.c.) + UY^ n^^n.i 

icr i^j^a i 

-ahjo'^ni^{ai + al) + '^hajo (afui + ■^j . (2) 

The meaning of the various operators and coupling con- 
stants is the same as in the Hamiltonian (1), except for 
the value of a which is equal to ^a. 

The electron-phonon interaction with coupling 
strength a is taken to be local, in view of the interaction 
with primarily local phonon modes which favours small 
polaron formation. U denotes some effective Coulomb 
repulsion - not to be confused with a Hubbard-type on- 
site interaction - on the effective sites, which, we stress 
again, ought to be considered as being composed of small 
molecular clusters such as a diatomic molecular unit or 
similar. Electrons described by this Hamiltonian can 
gain energy either by becoming itinerant and remaining 
essentially uncoupled to the lattice or, alternatively, 
getting self-trapped. In the first case the energy gain is 



described by the intrinsic band dispersion 

ek = i?(l-7k) j^^-Y,e^\ (3) 

^ s 

6 denoting lattice vectors linking nearest-neighbor sites. 
For self-trapped polaronic states, the local part of the 
Hamiltonian (^ can be transformed by a shift in the 
phonon variables a\'^^ a|^^ +Q:(n,;| +nii) upon which 
it acquires the form 

Hioc (£> - ^ - £p) riia- - (2£p - t/) nj^riii 

ia i 

+ Y^hwo(a+a, + ^j , (4) 

with the polaron binding energy Sp = a^fiwo denoting 
the gain in energy by the corresponding polaron level 
shift. The effective on-site interaction given by (2ep — [/) 
controls the possible formation of bipolaronic states 

p+X^^4c+e--(--?) (5) 

when 2£p — C/ > 0, with the single polaron states lying 
above such bipolaron states. For 2ep — U < we have the 
inverse situation with the bipolaronic level lying above 
that of the single polarons. What we are interested here 
is the first case. 

Let us begin by considering the situation of sufficiently 
strong interaction such that the bipolaronic level lies be- 
low the bottom of the unrenormalized electron band, i.e., 
for 2£p — \V > D, as illustrated in Fig. 1. In that case 
bipolarons can acquire itinerancy by decaying into vir- 
tual states of pairs of itinerant electrons and reassem- 
bling subsequently on some neighbouring site. Within 
the usual Lang-Firsov approximation, neglecting phonon 
creation as well as annihilation processes during the 
charge transfer from one site to the other, the correspond- 
ing effective bipolaron hopping integral is then given by 
t** ~ t2g-2a=/(2£p - iJ7 - !))§. Within this approxi- 
mation one obtains a system of itinerant bipolarons on 
a lattice which, at low temperatures, can condense into 
a superfluid state, giving rise to what has been termed 
bipolaronic superconductivity. This theoretical proposi- 
tion has ever since remained a matter of dispute, ques- 
tioning if such a normal-state itinerant Bose liquid can be 
achieved in real materials. The main arguments against 
it is that the stability of local bipolarons requires a rel- 
atively high value of > {D + ^U)/2hwo for realistic 
values of D, U and ojq and hence leads to exceedingly 
small values of the bipolaron bandwidth. A further ham- 
pering factor for the realization of such bipolaronic itin- 
erant band states is related to the retardation effects be- 
tween the motion of the charge carriers and the lattice 
deformations which surround them. This gives rise to 
dephasing between the two, as shown by exact diagonal- 
ization approaches Q and DMRG studies [Q, indicating 
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that incoherent processes in the optical conductivity are 
dominant. All these results suggest that coherent bipo- 
laron motion is most unlikely to occur in real systems and 
it was for that reason that the Boson- Fermion model was 
initially proposed]^. The aim was to capture the physics 
outside the regime of the extreme strong coupling antia- 
diabatic limit and to describe the intermediate coupling 
case, where the bipolaron level lies inside the band of the 
itinerant electrons (see Fig. 2) and where the supercon- 
ducting state is controlled by the phase correlations of 
these resonant bosonic bipolaron states. 

Quite generally, this model can be considered as a 
paradigm for cross-over situations. For the specific case 
we are interested here, it describes a cross-over between 
a BCS-type superconductor of Cooper pairs in the weak 
coupling regime and the hypothetical, yet to be exper- 
imentally verified, bipolaronic superconductivity in the 
strong coupling regime. In a totally different physical sit- 
uation, namely that of a strongly correlated electron sys- 
tem such as described by the Hubbard model, it has been 
recently shown that the intermediate coupling case can 
be mapped onto such a Boson-Fermion scenario §. Such 
an effective Boson-Fermion Hamiltonian was also derived 
recently for the exchange interaction between spinon sin- 
glets of the RVB electron-pairs and holonsQ]. The ques- 
tion to what extent such a BFM can emerge from a 
quite general class of interacting fermion systems, has 
been addressed from the point of view of a bosonization 
procedure of an intrinsically fermionic system]^. The 
presently much discussed molecular Bose-Einstein con- 
densate, involving entangled atoms in squeezed states has 
also recently been analyzed in the framework of such 
a BFM. 

A derivation of the BFM for moderately strongly cou- 
pled electron-phonon systems, starting from the Hol- 
stein model, is so far not available. In order to pro- 
vide the physically intuitive picture which initially had 
led us to propose this model, let us consider the case 
where the local bipolaron states are stable relative to 
single polarons. This eventuality can only be decided 
on the basis of direct experimental evidence for such lo- 
cal bipolaronic entities. The classical cases which ini- 
tially had led to consider such a situation are systems like 
amorphous semiconductors , Ti4O7 0, Naa;V205 [p^ , 
W03-3;fP^ and double valence fluctuation systems] 17 



like PbTe(Tl), just to name a few of them, and all of 
which represent such bipolarons as exclusively localized 
self-trapped states. 

When these systems are doped, the electrons respon- 
sible for pairing are at the same time responsible for the 
- albeit diffusive - transport. It is hence tempting to hy- 
pothesize that such bipolaronic states can occur as reso- 
nant states inside a band of itinerant electrons. In more 
complex systems, such as ternary and quaternary com- 
pounds, different electrons come into play. Some of them 
are prone to form bipolarons and others will remain itin- 



w 




FIG. 1: Schematic plot of the case of the bipolaron level (BP) 
falling below the band of itinerant electrons. 
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FIG. 2: Schematic plot of the case of the bipolaron level (BP) 
falling inside the band of itinerant electrons. 



erant, but being coupled together via some hybridization 
term. In systems like that, such a Boson-Fermion sce- 
nario imposes itself more directly. It is the eventuality of 
bipolaronic resonant states rather than localized states, 
which was at the origin of our proposal of the Boson- 
Fermion scenario iQ. 

The level scheme for the local bipolaronic states inside 
the band of the bare itinerant electrons for this situation 
is depicted in Fig. 2. Filling up such a system gradually 
with electrons, it is evident from Fig. 2 that, as long as 
the Fermi level remains below the level of the local bipo- 
laronic states, the charge carriers will be essentially bare 
itinerant electrons. But upon further increasing the con- 
centration of charge carriers, the itinerant electrons will 
eventually overlap in energy with the intrinsically local- 
ized bipolaronic level. 

The major supposition leading to this BFM is now to 
assume that this overlap can be simulated by some effec- 
tive on-site hybridization describing a charge exchange 
between the bipolarons and pairs of itinerant electrons 
in the form of 

vY.[ptX+a,^c,^- + p-X-c+clH , (6) 

i 

with the itinerant electrons and the locahzed bipolarons, 
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Eq. (5), treated as commuting fields. Such a supposition 
is supported by the fact that the spectrum of electrons, 
moderately strongly coupled to local modes, is known 
to be given by the bare spectrum for all k-vectors, ex- 
cept where their energy crosses the level of localized self- 
trapped electrons. In this regime of energies the elec- 
tron spectrum is strongly modified, as shown by ex- 
act diagonalization studies |l^ as well as variational 
methods [p^, with the spectrum becoming rather flat and 
approaching the local level of the self-trapped polaronic 



states. These are features which are characteristic of sys- 
tems where itinerant charge carriers are hybridized with 
localized charge carriers, and this was the center of our 
initial assumption of an effective hybridization between 
localized bipolarons and pairs of itinerant electrons, as 
given by Eq. (6). 

On a local level the dynamics of a polaronic system, 
involving two-particle resonant states, is controlled by 
the atomic limit of such a system, which we thus propose 
to describe by the following Hamiltonian 



(-4£p + U -2ii) 



E + ^) + ^ E [Pt^tcric^^ + P^ X, C+C+] -I- E fi^O 



(7) 



iU - 2/i) E [pI + + " E [^^c,|C,| -t- c+ c+] - hw^a E (yPl + {a^ + + E ("^ + ^) ' 



with a bare bipolaron level given by As ~ U. 

The equivalent formulation of this Hamiltonian [the 
second part of Eq. (7)] is obtained after employing an in- 
verse Lang-Firsov-type transformation Hat U^HatU 
with U = e°'^°-'-°'t'^^P'+i\ Adding to this local Hamil- 
tonian the term describing the kinetic energy of the itin- 
erant electrons, leads immediately to the BFM Hamilto- 
nian given by Eq. (1). 



COMPETITION BETWEEN LOCAL AND 
GLOBAL COHERENCE 

THE PSEUDOGAP PHASE 

The basic physics contained in this BFM model orig- 
inates from a competition of the local electron pairing, 
acting as a prerequisite for superconductivity, with the 
delocalizing effect of the itinerant electrons which fa- 
vors non-local Cooper pairing as the temperature is low- 
ered and the superconducting phase is approached. The 
change-over from local pairing to Cooper pairing is not 
only manifest in the electronic properties of such systems 
but is expected to be accompanied by corresponding lat- 
tice deformations, the study of which is the issue of the 
present paper. 

We shall examine this feature for the two following 
extreme limits: 

i) the normal-state phase, characterized by the exis- 
tence of a pseudogap in the density of states (DOS) of 
the itinerant electrons, which will be treated here in an 
approximate fashion within the atomic limit of the BFM, 

ii) the superconducting phase, which will be treated on 
a mean-field level, taking into account the phase locking 



of the bound electron pairs. 

Throughout this study we choose as a representative 
example a set of parameters presenting a situation close 
to the fully symmetric case. This implies an effective 
bipolaron level (corresponding to the energy level for the 
bound pairs) given by As — 4ep ~ 0, resulting in a fi- 
nal state bipolarons level lying close to the center of the 
band of the itinerant electrons, after the boson-phonon 
coupling has been switched on. We moreover want to 
consider the case ritot — 2np + ub = 2, giving a con- 
centration of electrons np ~ 0.5 and a concentration of 
bound pairs — 0.5 and choose v = 0.25, coq = 0.1 and 
a = 2. 

Let us start from the purely local physics, described 
by the atomic limit {t = 0) of the BFM corresponding 
to the Hamiltonian Hat, Eq. (7). This limit has been 
studied by us previously in connection with the pseu- 
dogap phenomenon and the incoherent background in 
the electron spectral function, which, in such a scenario, 
arise from phonon shake-off processes. For the choice 
of parameters adopted here, the physics controlling the 
change-over from the normal into the superconducting 
phase can essentially be represented by the one-, two- 
and three-particle eigenstates of Hat 

\l,irj = c+|0)|0)|Z> (8) 

I2,;r* = E(<44+^?n'^^)|o)|o)i"> (9) 

n 

with |«, j}"* denoting the z-particlc j-th eigenstate. The 
opening of the pseudogap at a certain temperature T* 
shows up in the rapid spectral weight increase, with de- 
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FIG. 3: Evolution of the fermionic DOS with temperature, 
showing the opening of the pseudogap at the chemical poten- 
tial (w = 0), as the temperature decreases below T*. 



creasing temperature below T*, of the scattering ampli- 
tude involving transitions from the lowest energy two- 
particle bonding state |2, 0}°* to the lowest energy single- 
particle states 1 1,0};^' and from the lowest energy three- 
particle states |3, 0}J^* to |2,0}°*. These contributions to 
the DOS of the itinerant electrons are manifest by the 
wings (see Fig. 3) for energies below and above that of 
the single particle non-bonding contributions, which have 
spectral weight around the chemical potential (equal to 
w = 0) and which arise from transitions not involving any 
of the two-particle states. Decreasing the temperature 
below T*, the spectral weight at w = decreases rapidly 
and transfers to the wings of the DOS, thus resulting 
in the opening of a pseudogap. Taking into account the 
dispersion of the itinerant electrons in an approximate 
fashion, such as within a Coherent Potential Approxima- 
tion (CPA) approach [p^, the main features seen in Fig. 3 
are maintained, simply leading to a smearing of the DOS. 
The two-particle eigenstates of Hat, Eq. (9), are obtained 
by expanding the local oscillator states in terms of the 
eigenstates \n >i of the undeformed harmonic oscillator 
at a given site i. |0) and |0) denote respectively the vac- 
uum states of the itinerant electrons and of the bound 
electron-pairs on such a site. The coefficients and v^^ 
are determined by exact diagonalization of this atomic 
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FIG. 4: Comparison of P{n) in the atomic limit and in the 
superconducting state for different temperatures T. 



limit problem within a restricted Hilbert space of those 
excited oscillator states. It is generally sufficient to take 
into account up to 50 such states in order to get converg- 
ing results. 

The qualitative changes in the lattice deformations ex- 
pected to occur when going from the pseudogap phase 
into the superconducting phase are essentially described 
by the lowest energy one-, two- and three-particle states, 
which contain the information on the relative weight of 
the bound electron-pairs (vf*) and of the induced pairing 
in the subsystem of the itinerant electrons (u°^). The 
first ones are intrinsically coupled to the local lattice de- 
formations while the second ones acquire such local lat- 
tice deformations due to the charge exchange mechanism 
acting between the two. As we shall see below, the ef- 
ficiency of this transfer of polaronic features is notice- 
ably different in the normal and in the superconducting 
phase. It shows up in the composition P{n) of the co- 
herent phonon states describing the local phonon clouds 
surrounding the charge carriers on a given site i, which 
in the normal state (atomic limit) is given by 



Pat{n) = ""'{mJla+almjr' 



(11) 



Z denoting the partition function and E'^^^i the eigen- 
values of Hat associated with the eigenstates jm, Z}"*. 
Upon decreasing the temperature, Pat{n) shows a re- 
distribution of weight from high to small values of the 
phonon number n [see Fig. 4(a)]. In the temperature 
regime of interest here, i.e., from slightly above T* all 
the way down to T = 0, the contributions to P{n) come 
essentially from the lowest energy one-, two- and three- 
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FIG. 5: Decomposition, for different temperatures, of P{n) 
in the normal state (full squares) into contributions coming 
from the one-, two- and three-particle states (Pi(n), P2{n) 
and Psin), respectively). Piin) and P2{n) represent the con- 
tributions coming respectively from the lowest energy two- 
particle state (bonding) and from the first excited two-particle 
state (antibonding). 

particle states, Eqs. (8-10). As illustrated in Fig. 5, these 
various contributions have noticeably different features 
as far as the phonon distributions are concerned. The 
one-particle states are totally decoupled from the lattice 
and hence have a contribution only for n — Q, while the 
three-particle contribution are maximally coupled to the 
lattice showing a broad peak in this distribution which 
corresponds to that of a shifted oscillator ground state. 
The two-particle bonding state |2,0}'^* has features of 
a displaced oscillator which is particularly noticeable at 
low temperatures where it is the main contribution, to- 
gether with a much reduced contribution from the two- 
particle antibonding state |2, 1}"*. That latter contri- 
bution to P{n), considered without thermal weighting 
factor, clearly shows strong contributions for low values 
of n, which is indicative of weak correlations between the 
electronic degrees of freedom and those of the lattice. 
The onset of the local coherence between bound pairs 
and of the induced pairing in the itinerant subsystem 
can be tracked in the pair distribution function (PDF) 
g{X), given in the normal state (atomic limit) by 

gat{X) = 1 ^ e-''^- f {m, 1\6{X - X,)|m, Ijf , 

(12) 

which describes the probability of finding a local defor- 
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mation of size X and which can be studied by EXAFS 
and pulsed neutron scattering experiments [pO| . gatiX) is 
obtained by expressing the eigenstates of Hat in terms of 
a real-space representation involving the harmonic oscil- 
lator wavefunctions. In Fig. 6(a) we illustrate the tem- 
perature evolution of the PDF. Its decomposition into the 
contributions arising from the one-, two- and three- par- 
ticle states given in Eqs. (8-10) are illustrated in Fig. 7. 
The one- and three-particle states contribute to the PDF 
in form of peaks centered around X = and around 
some finite value of X, respectively. This expresses the 
fact that the one-particle states (8) are completely de- 
coupled from the lattice, while the three-particle states 
(10) contain a bound electron-pair and thus give a max- 
imal shift in the PDF. The two-particle bonding and 
anti-bonding states, Eq. (9) with / ~ 0, 1, respectively, 
give contributions to the PDF with two peaks lying very 
close together, which is a result of the strong correla- 
tion between the bound electron pairs and the induced 
pairs in the itinerant electron system, characterizing the 
normal-state phase. At high temperatures, the PDF is 
thus given by a very broadened smeared out double-peak 
structure. But as the temperature is lowered below T*, 
local coherence between the bound electron pairs and the 
pairs of the itinerant electrons sets in. As a result, the 
opening of the pseudogap in the local density of states 
is accompanied by a change-over in the PDF from this 
very smeared out double-peak structure to a rather well 
defined single-peak structure. The latter is essentially 
due to the lowest energy bonding two-particle eigenstate 
contributing to a PDF with a maximum at some finite 
value of X. This is a signature of the strong correlations, 
which gradually build up as we go below T* , between the 
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bound pairs and the induced pairing in the itinerant elec- 
tron subsystem leading to roughly equal weight coming 
from the coefficients uf^ and vfjj^. 

THE SUPERCONDUCTING PHASE 

Let us now put such single units, as described by Hat, 
together in an infinite lattice, and let them interact with 
each other via the hopping term of the itinerant electrons. 
When the system undergoes a transition into a superfluid 
state of the charge carriers, the on-site correlations be- 
tween bound pairs and pairs of itinerant electrons will 
be weakened, since the k-space pairing of the itinerant 



electrons enters into competition with the strong on-site 
pairing which controls the physics in the normal state. 
As a result, upon entering the superconducting phase, a 
substantial modification of the coherent phonon states of 
the individual units will occur. 

We shall, for a moment, make the assumption that 
the deformation of the medium surrounding the bosonic 
charge carriers is not influenced by boson-fermion charge 
fluctuations. In that case the coupling of the bound 
electron-pairs to the lattice deformations can be elimi- 
nated, along the fines usually adopted in studies of the 
polaron problem, leaving us with an effective Hamil- 
tonian given by the first four terms of the Hamilto- 
nian (1) and with a renormalized hybridization coupling 
ve~°' Such a Hamiltonian has been studied in de- 
tail within different schemes, such as self consistent dia- 
gramatic techniques dynamical mean-field theory |22| 
and renormalization group procedures ||2^ , and gives rise 
to the following physics. As the temperature is lowered 
below a certain characteristic temperature T*, pairing of 
the itinerant electrons sets in, leading to an emptying out 
of the single electron spectrum near the chemical poten- 
tial and thus to the opening of a pseudogap. Upon further 
decreasing the temperature these electron-pairs become 
itinerant p^, eventually condensing into a macroscopic 
superfiuid ground state of the form 

n["»e-**'/2-f «,e+'*'/V+]|0). (13) 

i 

{ui,Vi being some model dependent parameters) and ex- 
hibiting soundwave-like collective excitations as a conse- 
quence of the rigidity of the local phases (f>i . 

In order to study how such a superfluid state affects 
the surrounding lattice deformations, we now investigate 
this problem on the basis of a mean-field treatment of 
the full Hamiltonian with the coupling of the phonons to 
the bound pairs explicitly taken into account. To this 
purpose, we introduce the following order parameters for 
the charge carriers [p5| 

- = ^E(44)' p - ^^(pt+pi) (14) 

i i 

which we assume to be homogeneous in space. The cor- 
responding mean field Hamiltonian is given by: 



rr TT TT ^0 

Hmfk = Hp + Hb - vpx + 



Hb = (Ab-2m)^ 



2 



^[pf +P^]- ^^oa ^ [pI + i ) (a, + a+) + huo ^ a. 



(15) 
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On the level of this approximation, the effect of the 
charge fluctuations between bound pairs and pairs of 
itinerant electrons is diminished, being taken into ac- 
count in a global, spatially independent, fashion. The 
itinerancy of the intrinsically localized bound electron 
pairs, induced by the exchange coupling ^ is not 
contained in the present approximation and therefore we 
cannot account for collective soundwave-like excitations 
and the spatially correlated lattice deformations which 
result from them. Nevertheless, the present study does 
permit to investigate local excitations of the lattice, con- 
trolled by the local fluctuations of the order parameters 
to which they are coupled, and which can be related to 
a number of experiments. 

Considering itinerant bosons on a deformable lattice 
with weak boson-phonon coupling [^, leads to a hy- 
bridization of the Bogoliubov modes (coming from the re- 
pulsive interaction between the bosons) and the Einstein- 
like phonon modes. At long wavelengths this results in a 
modified sound velocity of the Bogoliubov mode, with lit- 
tle spectral weight carried by the phonons, and a moder- 
ately renormalized Einstein mode, principally of phonon 
character. In such a context, what the present mean- 
field study permits us to do, is to investigate the effect of 
the superfluid state on such an Einstein-like mode. Due 
to the strong coupling between the bound electron pairs 
and the phonons, such an Einstein mode will undergo a 
substantial renormalization, as it will be shown below. 

The fermionic part Hp of TJmfa, being dynamically 
decoupled from the bosonic part Hb and hence from 
the phonons, can be cast into the standard form Hp = 
X]k^k(/3) (ttfcCkfc + PkPk)- The a's and /3's denote the 
operators of Bogoliubov quasi-particles having the dis- 
persion ek(p) = ± -y^ik^^/IP+T^^pIV^ which differs 
from the bare electron dispersion £k by showing a gap of 
size vp. The corresponding eigenstates of Hp are given 
by the standard BCS expression 



BCS\ 



n[ 



+ 



fcT -fci 



|0) 



(16) 



The bosonic part Hb oi Hmfa, showing intrinsic coupling 
between the bosonic charge carriers and the phonons, is 
diagonalized in terms of a set of states of the form 

I ar = E K(') + pt] io).i^ >^ (17) 



having eigenvalues Ei^{x). The eigenstates of the mean- 
field Hamiltonian, Eq. (15), are thus given by 



(18) 



As in the atomic limit, the diagonalization of Hb is 
carried out in a truncated Hilbert space spanned by the 
eigenstates |n > of the undeformed harmonic oscillator. 
We ultimately arrive at the following set of self-consistent 



equations determining the temperature dependence of 
the order parameters 

/3i?k(p) 



vp \ - 1 , 

TT7 / , , , tanh 
4^ ^ £k(p) 2 



<<exp [-/3i?r(a^)] , (19) 



In 



Z — '^iGxpl— PEi^{x)] denoting the partition function 
corresponding to Hb- This set of equations must be 
solved with the constraint that the total density of bosons 
and fermions remains fixed, i.e. ritot = np + 2nB 
where 



np 



-E( 



N 

1 
2 

1 
2 



'kcrCko-; — 



tanh^ 



Skip) 



i 



(20) 



In 



Let us now examine the modification of the composi- 
tion of the local coherent phonon states on the individual 
lattice sites (discussed above on the basis of the atomic 
limit for the normal state) when a macroscopic phase- 
coherent state of the charge carriers emerges and the 
bound elelctron-pairs condense into a superfluid state. 
The probability distribution of the phonons in the super- 
conducting phase is then described by 



I 



(21) 



The evaluation of Psc{n) in the superfluid phase fol- 
lows closely the calculation for Pat{n) in the atomic 
limit. In fact, the secular equation determining the co- 



efficients 



and 



'^In 



is formally the same as in the 



atomic limit, with the states c^'c^' |0 > (8) |0) (8) |« > 
and |0 > (g)p^|0) ®\n> spanning the Hilbert space be- 
ing now replaced by |0) ® \n > and p'''|0) ® \n >, re- 
spectively. The only modifications which then appear 
in this set of equations are (i) a renormalized hybridiza- 
tion constant v = vx, instead of v in the atomic limit, 
and (a) an energy level for the diagonal elements involv- 
ing the states |0) ^ \n > being equal to fiujon, instead 
of 2(1? — /i) -t- hujQn for the diagonal elements involving 
the states c^'c^ |0 > (8) |0) (8) |n >. This renormalization 
of V is important since it reduces the bare v by roughly 
an order of magnitude (for the set of parameters cho- 
sen throughout this paper), expressing the fact that the 
superfluidity leads to a considerable reduction of the on- 
site correlation between the configurations with bound 
electron-pairs present and with such pairs being absent. 
This has significant consequences on the composition of 
the local phonon clouds surrounding the individual sites 
and thus on their corresponding PDFs. 
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Just below the superconducting Tc (~ 0.0067 in our 
present case), Psc{n) is determined by an incoherent su- 
perposition of essentially the two lowest eigenstates of 
the mean- field Hamiltonian, Eq. (15), one giving rise to 
a broad peak at some finite n and the other to a sharp 
maximum at n = 0. At low temperatures, however, the 
ground state alone is sufficient to describe Psc{n). The 
contribution to Psc{n) leading to a peak at some finite 
value of n (associated with the coefficients Vq"^) arises 
from the deformation induced by the presence of the 
bound electron pairs, while the contribution to Psc{n) 
leading to a maximum for n = (associated with the 
coefficients wg^) arises from the absence of such an in- 
duced deformation for configurations where no bound 
pairs are present, with the oscillators being essentially 
undeformed. The coherent superposition of deformed 
and undeformed oscillator states in the ground state is 
a direct consequence of the macroscopic phase locking 
of the composite bosonic states in the superconducting 
phase, resulting from the variational ground state of Hb 
of the form 

Y[ i^e-'^^/^Xr + e^^'^X+p+') |0), ® |0 >, . (22) 

i 

The operators Xf^ are defined as 

^rio>^ = Y.<Mn>^ (23) 

n 

Xt\Q>^ - Y.<n{i)\n>^ (24) 

n 

and differ significantly from the usual polaron shift opera- 
tors = gT"(ai-ai )^ characterizing the deformations of 
a simply shifted oscillator state. Macroscopic phase lock- 
ing and homogeneity imply that the coefficients UQ'^{i) 
and t'ora(*) well as phases are the same for all sites. 
Eq. (22) is a direct generalization of the superfluid ground 
state of a system of hard-core bosons on a lattice, which 
is of the form given in Eq. (13), and thus shows explicitly 
the locking together of the charge {pf) and lattice (Xf) 
degrees of freedom on a local level. We shall next address 
ourselves to the question of how such a macroscopic co- 
herent quantum state of the phonons can be tracked ex- 
perimentally and, in particular, what difference there is 
to be expected with respect to the normal state features 
of the PDF, examined in the previous section. The PDF 
in the superconducting phase is given by 



T=0.001 



gscix) 



-f3Ei''(x) sc 



mx - xmi 



(25) 



In Fig. ^ we compare the PDF in the normal state and 
in the superconducting phase for different temperatures. 
The features examined above on the basis of the compo- 
sition of Psc{n) are perfectly refiected in the PDF. In the 
normal state one observes a single peak at some finite 



superconducting phase P(n) □ 

T,=0.00665 p„(„) 

Pi(n) - -o - 



0^ 




FIG. 8: Decomposition, for different temperatures, of P(n) 
(squares) in tfie superconducting pfiase in terms of tfie con- 
tributions Pq (n) coming from the ground state (triangles) and 
Pi(n) coming from the first excited state (circles) of Hb- 



Q 




FIG. 9: Decomposition, for different temperatures, of g{^) in 
the superconducting phase (full lines) into the contributions 
<7o(C) coming from the ground state (dashed) and 51(C) coming 
from the first excited state (pointed) of Hb- 



X, indicative of well defined shifted oscillator states. In 
the superconducting phase, on the contrary, one observes 
a weakly temperature dependent double-peak structure, 
which, as we can see from Fig.|^, arises essentially from 
two contributions associated with the ground state and 
the first excited state of Hb- At temperatures just below 
Tc this double-peak structure results from an incoher- 
ent superposition of these two states, corresponding to 
a situation where localized bipolarons are randomly dis- 
tributed over the lattice. The statistical distribution of 
empty and doubly occupied sites will thus give rise to 
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that double-peak structure. In the hmit of low tempera- 
ture, on the contrary, the PDF arises exclusively from the 
ground state, with a peak centered near X ^ (being 
determined by the coefficients ug^J and a peak centered 
around a finite value of X (being determined by the coef- 



with 



ficients 



The appearence of a double-peak structure 



in the PDF has long since been recognized as a signa- 
ture of delocalization of the charge carriers in polaronic 
systems] 27 . The present study shows that, as the super- 



conducting state sets in and is strengthened upon low- 
ering the temperature, the macroscopic phase coherence 
of the system not only involves the charge carriers but 
also the deformations which surround them. This is only 
possible if at T = these deformations are controlled by 
a single quantum state, i.e., the ground state JJ - \0}l''. 
The superconductivity-induced delocalization of the self- 
trapped entities which characterize the normal state, is 
caused by the locking together of the local lattice de- 
formations in a coherent superposition of deformed and 
essentially undeformed oscillator states, resulting in the 
developement of a double-peak structure in the PDF. At 
least in the scenario presented here, this result does not 
bear out any undressings of the charge carriers, strip- 
ping off partially the phonon clouds, as one enters the 
superconducting phase. Such an effect would show up 
in a PDF which, in the superconducting phase, would 
be given by essentially a single-peak structure centered 
around X ^ 0. 



DYNAMICAL LATTICE PROPERTIES 

The discussion in the previous Section demonstrated 
the qualitative difference which exists between some ba- 
sic lattice properties in the normal and the supercon- 
ducting phase, and largely being controlled by the in- 
duced electron pairing in the subsystem of the intrinsi- 
cally uncorrelated itinerant electrons. What results is 
a competition between the local pairing, characterizing 
the normal-state phase, and the k-space pairing in the 
superconducting phase, expected to be particularly rel- 
evant around the superconducting transition where the 
spatial phase fluctuations play a major role and residual 
local electron pairing |p9|l persists in the superconducting 
phase. The qualitative difference in the structure of the 
PDF in the two phases leads us to expect squeezing ef- 
fects of the local coherent phonon quantum states which 
are stronger in the pseudogap phase than in the super- 
conducting phase. Associated with it, one should expect 
larger fluctuations of the local deformations of the envi- 
ronment and of the momenta of the ions representing it. 
These fluctuations are given by 



Xr - 



P, 



2MuJo 
hM luq 



(a+ + a,) 



(27) 



The expectation values in Eq. (g^ have to be calcu- 
lated with respect to the eigenvectors of Hat and Hb 
for, respectively, for the normal and the superconduct- 
ing phase. The results for these fluctuations in the two 
phases are illustrated in Fig.|l^. In particular, the prod- 
uct AXfAP^ clearly shows that the coherent phonon 
states in the normal state are considerably squeezed, 
i.e., this product is for T — > very close to the lower 
limit of the Heisenberg uncertainty principle, equal to 
h?/A. In the superconducting phase, on the contrary, we 
observe practically no squeezing effect with decreasing 
temperature, the product AXfAP^ being always big- 
ger than what would be expected in the normal state in 
this temperature regime. A possible experimental ver- 
ification of these features comes from neutron absorp- 
tion measurements which measure the kinetic energy 
Ekin — ^{Pi) of the local lattice vibrations. Indeed, 
the experimental results shown that, as the temperature 
is reduced, Ekin saturates in the superconducting phase 
and is distinctly above the values expected in the nor- 
mal state at the same temperature. This fact is born out 
in our study, as can be seen by inspection of Fig.[lO|(b). 
Finally, an issue which ought to be of crucial interest in 
studying the lattice properties connected with the onset 
of the superconducting order or of the electron pairing 
at T* in the normal state, is the question of the changes 
in the phonon spectral properties upon going from one 
phase to the other. Strong changes in the frequency of 
specific local modes, induced by strong electron-phonon 
coupling, have been well established by photoinduced in- 
frared absorption measurements ||30| and been related to 
strong local displacements of the ions. In this context, 
the appearence of new additional modes occuring with 
changes from the normal into the superconducting phase 
have also been reported, although this is still a matter 
of debate ||3l|, More established are the strong in- 

tensity changes of certain Raman active modes and the 
strong frequency renormalizations and lifetime broaden- 
ing occuring at such transitions js^, Such effects can 
best be studied by energy loss spectra in inelastic neu- 
tron scattering or Raman spectroscopy where phonons 
are emitted. The cross-section for it is given by 



^mi(c^) = J die'"*^emi(i) 



with 



AX^ = (Xf) - (X, 



AP^ = (P^ 



(26) 



Aemi(i) a{0)a\t) > 9{-t) 



(29) 
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FIG. 10: Comparison of (a) the positional and (b) the mo- 
mentum fluctuations of the lattice in the normal (N) and in 
the superconducting (SC) phase as a function of relative tem- 
perature T/To (To = {TcT*}) and with = {l/2)MLj^£p. 
(c) Degree of the squeezing of the coherent phonon states as 
a function of T/To- 



Here ei = {Ef,E^''{x)} and \l >= {\mj >'^M? >^'=} 
are the eigenvalues and eigenvectors of the Hamiltonians 
Hat and Hb, respectively. We illustrate the results of 
such phonon emission spectra for the two phases in Figs. 
|lT| and 12. Quite generally, one observes a quasi-elastic 
contribution at frequency zero, which is the standard sig- 
nature of coherent phonon states. As the temperature 
is lowered below either T* or Tc and the itinerant elec- 
trons get correlated by local pairing in the normal state 
or by Cooper pairing in the superconducting phase, the 
local phonon mode first splits into several modes and fi- 
nally evolves into a well defined mode which is softer in 
the pseudogap phase and harder in the superconducting 
phase. These features of the phonon spectral properties 
are expected to hold qualitatively true. A quantitative 
link-up between the normal phase and the superconduct- 
ing phase as the temperature is decreased, would require 
a full description of the local boson-fermion exchange 
correlations, incorporating the electron itinerancy in the 
normal state and allowing for local electron pairing in 
the superconducting phase. These are presently largely 
unresolved fundamental problems in such system where 
amplitude and phase order occur independently, as in 
the high-temperature superconductors with its pseudo- 
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FIG. 11: Evolution of the phonon spectrum in emission ex- 
periments as a function of temperature in the normal state. 
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FIG. 12: Evolution of the phonon spectrum in emission exper- 
iments as a function of temperature in the superconducting 
phase. 



gap phase above Tc where local electron pairing persists 
in an incoherent fashion. 



CONCLUSIONS 

The present study is an attempt to examine the con- 
densed state of a system of composite particles such as 
charge carriers surrounded by clouds of bosonic excita- 
tions of the neighboring environment to which they are 
coupled and to see to what extent a macroscopic coher- 
ent quantum state can be induced in the system of these 
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bosonic excitations. As a specific example we chose a 
system of electrons, moderately strongly coupled to local 
lattice deformations, which can be described in terms of 
resonant local-pair states (localized bipolarons) inside a 
Fermi sea of essentially uncorrelated itinerant electrons 
with which they interact via a charge exchange term. We 
find noticeable qualitative differences in the structure of 
the local coherent quantum states of the phonons in the 
superconducting phase and in the phase above it which 
corresponds to a normal state characterized by a pseu- 
dogap in the DOS of the single-particle spectrum. Sev- 
eral measurable quantities related to static and dynamic 
lattice properties were examined, such as the pair distri- 
bution function, the fluctuations in the lattice positions 
and the corresponding ion kinetic energies of the atoms 
making up such a lattice and, finally, the phonon spec- 
tral properties given by energy-loss spectroscopy (mea- 
surable by inelastic neutron scattering or Raman spec- 
troscopy). The difference between local electron pairing 
in the normal state and k-space pairing in the supercon- 
ducting phase has noticeable consequences on such mea- 
surable quantities, resulting in (i) a single-peak struc- 
ture of the PDF in the normal state and a double-peak 
structure in the superconducting phase, (ii) distinctively 
larger positional fluctuations and kinetic energy of the 
atomic clusters (making up the lattice) in the supercon- 
ducting phase than in the normal state and, (iii) asso- 
ciated with it, a more effective squeezing of the phonon 
coherent quantum states in that latter. For the supercon- 
ducting phase, the present mean-field study is restricted 
to local phase fluctuations of the superconducting state 
given by a generalization of a phase- locked state of bosons 
on a lattice, Eqs. (17, 18, 22), in which these bosons cor- 
respond to composite particles in form of bipolarons [see 
Eq. (5)]. The examination of soundwave-like collective 
excitations of such a system is planned for some future 
work. There, the question arises if the local dynami- 
cal deformations do get spatially correlated upon enter- 
ing the superconducting state, as possibly suggested by 
ion channeling experiments js^. A possible approach to 
solve this problem could be to integrate out the fermionic 
sector of the BFM Hamiltonian Eq. (1), resulting in a 
Gross-Pitaevskii Lagrangian|^ for the effective itinerant 
bosonic charge carriers coupled to the phonons, followed 
by a treatment similar to that adopted for bosons on a 
defromable lattice 
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